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ABSTRACT 



We construct a simple model for stationary, axisymmetric black-hole magneto- 
spheres, in which the poloidal magnetic field is generated by a toroidal electric cur- 
rent in a thin disk with the inner edge, by solving the vacuum Maxwell equations in 
q Schwarzschild background. In this work, to obtain a concise analytical form of the mag- 

netic stream function, we use the approximation that the inner edge is far distant from 
the event horizon. The global magnetospheric structure with the closed-loop and open 
field lines threading the inner and outer parts of the disk is explicitly shown, claiming 
that the model is useful as a starting point to study astrophysical problems involving 
inward disk-driven winds to a black hole and outward ones to infinity. The asymptotic 
shape of the field lines at the event horizon becomes nearly cylindrical, while at infin- 
ity it becomes conical. The magnetic spot in the disk connected with the black hole 
through the loop field lines occupies a very narrow region with the ring area roughly 
equal to the horizon area. By taking account of the existence of a uniform (external) 
magnetic field, we also obtain the model for collimated open field lines. Then, it is 
found that the magnetic connection between the black hole and the disk breaks down 
if the uniform field is strong enough. Considering slow rotation of the magnetosphere 
and angular momentum transfer by inward winds from the disk, the final discussion is 
devoted to gradual disruption of the closed loops due to radial accretion of disk plasma 
toward the black hole. 



Subject headings: black hole physics — accretion disks — magnetic fields 
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1. INTRODUCTION 

It is widely believed that a super massive black hole surrounded by an accretion disk can work 
as a central engine of active galactic nuclei. If magnetic field exists around the black hole and 
the disk, it is important to understand the magnetospheric structure for explaining the highly 
energetic phenomena. The magnetohydrodynamical (MHD) approach in general relativity may 
be necessary to study various features of the black hole magnetospheres. Unfortunately, even if 
the magnetosphere is assumed to be stationary and axisymmetric, any global models describing 
both inward winds in the strong-gravity region near black hole and outward winds in the far distant 
region have not been constructed, in spite of a large number of works devoted to the MHD problems. 
This is mainly because it is very difficult to analyze the highly nonlinear Grad-Shafranov equation 
for the magnetic stream function which defines the poloidal field lines (see e.g., Beskin (1997) for 
review) . 

To avoid the mathematical difficulty in the full MHD system, one may consider vacuum solu- 
tions of the magnetic stream function as a practical approximation applicable to the magnetically 
dominated models. In particular, the well-known examples such as split-monopole, uniform and 
paraboloidal fields in Schwarzschild geometry (Michel 1973; Wald 1974; Blandford & Znajek 1977) 
have been used as the zeroth-order approximation for the perturbation method in slow-rotation and 
force-free fields, from which the first-order perturbation has been derived for giving a weak toroidal 
magnetic field and discussing the Blandford-Znajek process (Blandford & Znajek 1977). Another 
important MHD problem to be studied is fluid motion along poloidal field lines. For examples, the 
existence of accretion flows onto the black hole passing through some MHD critical points has been 
shown under a fixed shape of field lines (Takahashi et al. 1990). The vacuum solutions allowing a 
field configuration of astrophysical interest will become a useful tool in an attempt to discuss the 
trans-critical motion of outgoing and ingoing magnetized winds, except in the region where the 
effect of fluid inertia should dominate to change crucially the field structure. 

Of course, even for the vacuum fields, some current distribution should be assumed to exist in 
a restricted region of the magnetosphere. In this paper we focus on the magnetic field generated by 
a toroidal current distribution in a thin disk around the Schwarzschild black hole. The source of the 
above-mentioned fields with split-monopole and paraboloidal structures can be also disk currents. 
However, such models contain only open field lines and fail to describe a magnetic connection of the 
disk with the black hole. This will be partially due to the absence of the disk inner edge separated 
from the event horizon. The magnetic connection by closed field lines is expected to exist and 
have astrophysically important effects (Punsly 1991a; Gruzinov 1999; Blandford 2000; Li 2000c; 
Krolik 2000). Hence, our purpose here is to present explicitly a global magnetospheric model with 
both closed and open field lines threading a thin disk and connecting magnetically the inner and 
outer parts, respectively, with the black hole and the far distant region, which is relevant to the 
problems of inward and outward disk-driven winds. Though the vacuum model would miss some 
important features of more realistic magnetospheres, our modeling could be a preliminary step for 
understanding the various MHD processes in further detail. 
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To obtain general models based on a disk current distribution, it may be convenient to consider 
linear superposition of the magnetic field generated by a single ring current located at a fixed radius 
on the equatorial plane which has been found in Kerr geometry (Petterson 1975). Assuming the 
current distribution from an inner edge to an outer one, Li (2000a) has derived a formal expression of 
the magnetic stream function in terms of the infinite sum of the multipole fields, without revealing 
the global shape of the field lines. Because the each multipole field of which the structure has 
been discussed by Ghosh (2000) still fails to give closed field lines regular at the event horizon, the 
infinite sum is essential to understanding the global shape. 

In §2 we would like to develop a different approach to arrive at the infinite sum formula for the 
magnetic stream function. In §3 we succeed in reducing it to an analytical simple form by locating 
the disk inner edge at a radius much larger than the horizon radius and the outer edge at infinite 
radius. The global structure of closed and open field lines is shown, and it is found that only the 
very narrow region of the disk separated from the inner edge can be magnetically connected with 
the black hole, together with the result that the field lines become nearly cylindrical at the event 
horizon and conical at infinity. In §4 we also consider a modification of the global shape due to an 
addition of a uniform external field. Though collimated field lines appear at the far distant region, 
it is found that for a strong uniform field the magnetic connection between the black hole and the 
disk breaks down. Finally we discuss the MHD effect due to slow rotation of the magnetosphere 
(namely, angular momentum transfer by magnetically dominated inward winds from the disk to the 
black hole) , which could lead to radial accretion of disk plasma and gradual disruption of the closed 
field lines (Punsly & Coroniti 1989). In the following we use the geometrical units with G = c = 1. 



2. THE MAGNETIC STREAM FUNCTION 



Let us consider a stationary, axisymmetric magnetic field in Schwarzschild background of the 
metric 

ds 2 = - (l - dt 2 + (l - dr 2 + r 2 (d9 2 + sin 2 0# 2 ) , (1) 

where M is the mass of the black hole. The poloidal components of the magnetic field may be 
given by 

W x 



Bp ~ 2vrrsin0 ' (2) 
where is a unit vector in the azimuthal direction, and *(r, 6) is the so-called magnetic stream 
function defining the poloidal field lines by * = constant. If the vacuum field is assumed, the 
Maxwell equations reduce to the simple form 
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which allows the separable solutions 



* = R u (r)Qu(6) • 



(3) 
(4) 
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with the separation constant v 

For the analysis of these functions it is convenient to introduce the variables 

T 

w = 1 , x = cos 6 . (5) 

M w 

Then, from the regularity at the polar axis 9 = 0, ir we obtain the angular part written by 

Gu = [ Pu(x')dx' , (6) 

J\x\ 

where P v is the Legendre function of the first kind. Though in this paper we focus on the field 
with equatorial symmetry, the parameter v is not limited to a positive odd integer for allowing 
a discontinuous change of the radial component of B p on the upper and lower surfaces of the 
equatorial disk. 

To determine the radial part R u , the condition for regularity at the event horizon r = 2M and 
at infinity r — > oo becomes important. We use the Legendre function P u (w) to give the solution 
regular at w = 1 of the form 

/w 
P v {w')dw' , (7) 

where we assume to be v = (—1/2) + ik with arbitrary positive k, because the oscillatory real 
function P u (w) has the amplitude decreasing in proportion to \j\fw at large w, and the linear 
superposition 

poo 

tf(r,0) = / f(k)R u (r)e u (9)dk (8) 
J o 

will be able to describe a global magnetic field satisfying the boundary conditions. 

Now we calculate the coefficient f(k) in the linear superposition corresponding to the case of 
a disk current with the inner edge at w = w q > 1. Recall that in previous works the solution was 
written by an infinite sum of the multipole fields with v of natural numbers. Then, in the distant 
region w > wq the Legendre function P v (w) should be replaced by Q v {w) with a change of the 
coefficient due to the condition for continuity of *(r, 9) at w = wq. We can avoid these slightly 
complicated steps by using the integral form (8). Though the toroidal surface current in the thin 
disk is proportional to r~ 2 {d^ /d9) at 9 = tt/2, here we rather focus on the equality 

d fd^\ w + l d , , , 

-h(w) , (9) 



- 



dr\d9 J M dw 

which is evaluated at the upper surface of the equatorial plane. The remarkable point is that the 
function h(w) is given by 

/•oo 

h(w) = / g(k)P u (w)dk , (10) 
Jo 

where we have 

g(k) = P v (0)f(k) , (11) 
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and the inversion formula (Erdelyi et al. 1953) holds as follows, 

/oo 
P u (w)h(w)dw . (12) 

Then, if the function h(w) is specified according to the current distribution, we can derive the 
coefficient f(k) giving the magnetic stream function *(r, 0) valid in the whole region. 

Further, to obtain the explicit form of f(k) through the well-known integral formula involving 
two Legendre functions, let us assume in the disk region w > wq to be 

h(w) = A Q s (w) . (13) 

Though Aq may be an arbitrary real parameter, we limit s to be positive for assuring that d^/dO 
does not diverge at large w, where the Legendre function Q s (w) of the second kind falls off as w~ s ~ l . 
(If any different model of the current distribution becomes necessary, the linear superposition of 
Q s with various s in the range < s < oo is a possible procedure to construct the solution <!'. 
However, we would like to emphasize that this simple example for the current distribution can be 
useful for revealing interesting features of the magnetospheric structure, as will be shown in the 
following.) Because in the inner region 1 < w < u>o, where the field lines should vertically thread 
the equatorial plane, we have 

h(w) = , (14) 



it is easy to arrive at the result 



_ A)fctanh(Trfc) , , 

9{k) ~ fc 2 + [ s+ (l/ 2 )]2^' ^ 



where 

F v = K 2 -l) 



P^) d ^P±-Q s { W *) dMm) 



(16) 



dwo dwo 

Strictly speaking, even if the condition (14) is assumed, equation (9) only implies that d^/d6 is 
constant at the upper or lower surface of the equatorial plane in the inner range 1 < w < Wo. 
Hence, we must confirm the validity of equation (15), by calculating the value on the event horizon 
w = 1, which is given by 

— = ^ 2g(k)dk , (17) 

at the upper surface of the equatorial plane. In fact, the right-hand side of equation (17) is shown 
to vanish by virtue of the equality 

~ f BT¥Tm? PAm)dk (18) 

which is a result of the inversion formula (12). 

We have obtained the magnetic stream function of the integral form (8) in which the coefficient 
f{k) = g{k)/Pv(0) is given by equation (15). To calculate the integral with respect to k with the 
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application of the residue theorem, it will be useful to consider the integration contour along the 
real axis (—00 < k < 00) and the large semi-circle in the complex fc-plane, based on the relation 
between two Legendre functions written by 

P v (z) = -cot\i{itk)[Q v {z) - Q.^iz)] . (19) 

7T 

If w < wo we replace P 1/ {wq) in the integrand according to equation (19), and the large semi-circle 
is chosen to be in the lower (or upper) half of the complex k-plane for the integral of the term 
involving Q v (or Q_„_i). Then, by taking account of the contribution only from the poles present 
at k = ±i{s + (1/2)] and at k = ±i[2n + (3/2)] where P„(0) vanishes for n = 0, 1, 2, • • -, we obtain 
the result valid in the inner region 1 < w < wo as follows, 

00 

* = ^ n G 2n+ iiW(r)6 2n+1 (0) , (20) 

n=0 

where 

_ ^o(-l)"(4n + 3) (2n + l)!! 
(2n + s + 2)(2n + l-s) (2n)!! 



and 



G 2n +i = (w - 1) 



dQ s (wo) , ,dQ 2n +i(w ) 

,^ iW -sr " Qs(wo) d W0 



(22) 



In the outer region w > wq, however, we must use equation (19) for P u (w) [instead of P u (wq)] 
with the same choice of the large semi-circle, and the contribution from the poles is divided into 
the two parts 

* = *i + ^ 2 (23) 



with the forms written by 

*i = $^a n e2n+i(0) F 2n+1 R 2 n+i(r) -G 2n +i / P 2n +\{w')dw' 



(24) 



* 2 = 4^rfi-^V(l + ? )@ s m(r) (25) 



n=0 

and 

pa V 2 

where the function R u (r) for any is defined by 

^(r) = (w + l)Q v {w)- / Q u {w')dw' , (26) 



'UIO 

and T(z) is the gamma function. 

It is easy to check that the continuity of these expressions (20) and (23) for \I> at the boundary 
w = wq holds, because through the application of the expansion in terms of the Legendre poly- 
nomials to the function u(x) defined by u = P s (x) for x > and u = — P s (|x|) for x < we 
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obtain 

(-l) n (4n + 3) (2n+l)!!„ , , 1^/1-sX / s\ , , 
g o ( 2n + S + 2)(2„ + l- 8 ) W = 7? r (— ) r ( 1+ 2j" W ' (27> 

Further, equation (27) is useful to show that for a positive s the asymptotic behavior of \I> at 
w — > oo becomes 

/•oo 

* ~ -A / Q s (v/)dv/ x , (28) 
where the dependence on x is given by 

^ (-l)"(4n + 3) (2n + l)!! f 1 
<X) = ^(2n+l)(2n + 2) (2n)!! j^ l+ i{x)dx . (29) 

By substituting the value of s = into the equality (27), we can calculate the summation to 
be v = 1 — \x\. Hence, we conclude that this black hole magnetosphere based on a disk current 
model has the asymptotically conical shape of magnetic field lines. This asymptotic structure is a 
consequence of the disk currents extending to infinite distance. Though such a current distribution 
is unphysical, the existence of open field lines will be useful to discuss the problem of outgoing 
winds from the disk to infinity in our modeling. 



3. THE GLOBAL STRUCTURE OF FIELD LINES 

Though the asymptotic behavior of * at infinity is clear, the expression written by the infinite 
sum of the multipole fields is still complicated for understanding the features of the field line 
structure near the inner edge and the event horizon. Hence, let us assume that the inner edge is far 
distant from the event horizon. For example, the inner edge of the disk around the Schwarzschild 
black hole may be located at the innermost stable circular orbit r = 6M on the equatorial plane. 
Then, we have wq = 5, for which the approximation wo > 1 will be roughly allowed. 

In the following we focus on equation (20), because the multipole sum should arrive at the 
same result even if equation (23) is used. By virtue of the approximation wo 3> 1, we obtain the 
dominant behavior of the Legendre function such that Q s (wq) ~ u^ -1 . Except in the region near 
the event horizon where w is of order unity, we can calculate equation (20) under the assumption 
1 <C w < wq, which leads to the approximate relation 

iWw)Q w („o) * l2n + i; in + 3) {^) ■ (30) 

Of course, for the angular part involving the integral of P2 n +i(x) we cannot use such an approxi- 
mation. Hence, our key step is to rewrite the Legendre polynomial through the integral formula 

1 f n 

P 2n +i(cos8) = - / (cos + i sinO cos ip) 2n+1 d(f , (31) 
n Jo 
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which allows us to obtain 



Vt ~ A Qs{w )w ! dx' - I d(pK(y) 

J\x\ I 71 J0 



The variable y and the function K are defined by 

w 

y = — 

Wo 

and 

oo 



' + % COS if \J\ — (x 1 ) 2 



(_l)"(2n + l) (2n + l)!! y2n+1 



2n + 2 + s (2n + 2)!! 

n=0 v ' 

for which the infinite sum becomes possible to give 



(32) 

(33) 
(34) 



K = 



1 

y 



l 



VT+v 2 



-(s + l)y 



-s-1 



r\s-l 



v/i + (y') 2 



dy' 



(35) 



This expression of K(y) is also valid in the outer region w > wq, and it is easy to verify that 
at infinity where j/ > 1 the magnetic stream function given by equation (32) can keep the same 
asymptotic form as equation (28). 

Next, let us give the magnetic stream function near the event horizon. Because w is of order 
unity, the term with n = dominates in the multipole sum (20) under the approximation wq 3> 1, 
and we have 

A Q s (w ) 2 2 



-(i-x 2 )(w + iy , 



4(a + 2)w 

which shows the cylindrical magnetic field written by ^ = nB^r 2 sin 2 9 with the strength 

A Q s (wq) 



Bo 



AttM 2 (s + 2)wo ' 



(36) 



(37) 



This is in accordance with the result claimed for general multipole expansions by King, Lasota, & 
Kundt (1975). However, if wq is of order unity, the dominance of the cylindrical field in the sum 
(20) apparently breaks down even near the event horizon, and a more complicated structure of field 
lines will become possible there. 

We hereafter fix Bq to be positive in equation (37). Then, for s > 0, the magnetic stream 
function becomes negative at infinity, where we have 



- 1 



47r J B (Mu>o) 2 (1 - \x\ 



(38) 



Because this change of VP occurs for all 9, a closed field line of VP = connecting the equatorial 
plane with the polar axis should exist in the magnetosphere. This critical field line can thread the 
black hole along the polar axis \x\ = 1 and divide the magnetosphere into the inner region with 
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closed field lines and the outer region with open field lines. Note that if equation (32) is written 
by the sum of the multipole fields we obtain the r-dependence at x = as follows, 



* ~ B (Mw) 2 Y 2(5 + 2) 
uv ' ^ 2n + 2 + s 



n=0 



(2n + l)!! 
(2n + 2)!! 



1 2 



W 
Wq 



2n 



(39) 



which claims ^ to be positive everywhere on the inner equatorial plane w < wo between the disk 
and the black hole. Hence, the critical field line of * = should thread the disk at some point 
which we denote by w = w c (or r = r c ). For example, if we consider the model of s = 2 with K of 
the form 



K 



If 




(40) 



the numerical estimation of equation (32) gives w c ~ 1.9u>o. (The critical field line also reaches to 
the polar axis at w c± 1.8w;o, along which w is approximately constant.) 



To give the magnetic stream function approximately valid in the whole range from w ~ 1 to 



w 



wo, we modify equation (32) which reduces to ^ — > ttBq(Mw) 2 (1 



x 



in the limit w <C wq. 



We note that ^ given by equation (32) is smoothly matched to the cylindrical field (36) which is 
valid near the event horizon, only by multiplying the factor (w + l) 2 /w 2 which becomes equal to 
unity in the region w 3> 1. Hence, the modified form of the magnetic stream function should be 



* = 4(5 + 2)^0(^ + 1) 



\w /w)M 2 I dx' 

J\x\ 



dipK(y) 



(41) 



with K given by equation (35). 



We have arrived at the main result (41) available as a model for describing the global structure 
of the black hole magnetosphere. The numerical example of magnetic field lines for s = 2 and 
wq = 10 is drawn in Figures 1 and 2, which are useful to see the magnetic connection between 
the disk and the black hole in the inner region and the conical shape at large r, respectively. The 
magnetic spot in the disk connected with the black hole is found to be a very narrow region. This 
is because ~ BqM 2 near the event horizon, while \I/ ~ Bq(Mwq) 2 on the disk unless w ~ w c . If 
the width of the magnetic spot is denoted by Aw c (or Ar c ), at w = w c we obtain 



B n M 2 



(99_ 

\dw 



Aw r 



(42) 



which leads to Aw c ~ 1/wo -C 1 (i.e., Ar c <C M). Because the area 2itr c Ar c of the magnetic spot 
is just of order of the horizon area 16nM 2 , no significant amplification of the magnetic field occurs 
near the event horizon. Steady inflows of the disk plasma to the black hole may occur along the 
magnetic field lines, starting from the magnetic spot of the disk which is apart from the inner edge. 
Then, the inflows reach to the polar region at a distance roughly equal to Mwq and finally fall to 
the black hole along the cylindrical field lines. 
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Plasma outflows from the outer part of the disk may also emanate along the open field lines 
of < extending to infinity. It should be remarked, however, that the co-existence of two 
topologically distinct zones in the magnetosphere is not due to a general relativistic effect of the 
black hole, because such a global structure is also possible for an appropriate superposition of a 
dipole and a uniform field, giving closed and open field lines, respectively. The role of the black 
hole is rather to prohibit the presence of a dipole field at the event horizon. Hence, the closed field 
lines should be sustained by disk currents only, which is the point essential to our modeling. 

In astrophysical magnetospheres various MHD processes are expected to become important for 
producing highly energetic phenomena. For example, in the region very close to the event horizon, 
strong gravity forces the plasma to accrete in radial direction and to bend the cylindrical field lines 
(Punsly 1991b; Hirotani et al. 1992). The magnetic field (41) generated by disk currents provides 
a starting point to analyze such MHD effects of ambient and disk plasma on the magnetospheric 
structure, which will be discussed in the next section. 



4. DISCUSSION 

If the plasma inertia is taken into account even under the condition of magnetic domination, 
the asymptotic collimation from the conical shape of field lines should occur in logarithmic scales 
of the cylindrical radius (Chiueh, Li, & Begelman 1991; Tomimatsu 1994). Apparently such an 
inertial back-reaction of the plasma on the poloidal field does not work in the vacuum fields. Then, 
we present a model with collimated structure, by adding the uniform field (which is a solution of 
the vacuum Maxwell equations) 

^ u = TrB u r 2 sm 2 e (43) 

to given by equation (41). Let us denote the new magnetic stream function by \P'(= + ^f u ). 
Because ^ is negative in the far distant region, B u should be chosen to be negative for allowing 
a smooth collimation of the conical field lines. If the uniform field is weak, namely Bq + B u > 0, 
the value of remains positive at the event horizon to keep the magnetic connection between the 
disk and the black hole (Fig. 3). However, if Bq + B u becomes negative, the field lines threading 
the black hole can reach to infinity only (Fig. 4). 

Of course, this addition of a uniform field is just a mathematical procedure, and the strong 
collimation shown in Figures 3 and 4 is not a consequence of the physical confinement mechanisms 
by self-generated toroidal fields and/or by external pressure with a boundary at finite cylindrical 
radius. It remains quite uncertain whether the details of the field- line shape shown in these figures 
become approximately valid in magnetically dominated wind regions or not. Nevertheless, it is 
interesting to note that the vacuum models can successfully describe the transition of the structure 
of field lines threading the black hole (i.e., from magnetic connection with the disk to that with 
the remote load), which may occur in more realistic magnetospheres as an observable change of the 
astrophysical activity. 
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To discuss the magnetic activity by using the vacuum field ty r (or $) as a background field 
perturbed by magnetically dominated winds from the disk, let us consider the magnetosphere 
slowly rotating with the angular velocity of a magnetic field line and the hole's one (i-e., 
rsmOQp <C 1 and MVLh — a /AM <C 1). If a magnetic field line threading the black hole is closed in 
the disk at r = r d (see Figs. 1 and 3), the angular velocity ftp is nearly equal to the local Keplerian 
one, i.e., 

n F ~ n K = {M/rl) l/2 , (44) 

and we have r d ~ Mxwq, where x is a numerical factor in the range 1 < x < 1.9 for the s = 2 
model. Now the approximation of slow rotation is assured by virtue of the estimation such that 
r^F — 1/ v 7 P w o ^ 1- However, the inner edge in realistic accretion disks cannot be too distant 
from the event horizon (namely, wq cannot be too large), and from the relation 

n H /n F ~ Mn H ((3w ) 3 / 2 , (45) 

we expect to be Qh < tip as a typical case for slowly rotating black holes. Hence, the inward winds 
from the magnetic spot carry angular momentum from the disk to the black hole, and the disk 
plasma which remains in the magnetic spot without flowing into the magnetic tube will be torqued 
and fall inward. This means that the disk currents located in the inner region and sustaining the 
loop field lines become unstable for radial accretion, and the closed loop field disconnected from the 
disk must disappear owing to the infinite redshift effect at the event horizon (Punsly & Coroniti 
1989). 

In the magnetically-dominated wind region the toroidal magnetic field Bp is approximately 
given by 

Bp — —AirrjL/aw , (46) 

where a and w are the lapse function and the cylindrical radius, respectively, in Kerr geometry. 
By virtue of the approximation of slow rotation we obtain a = \J\ — (2M/r) and w = rsm9. 
Along a poloidal field line the particle flux 77 per unit flux tube and the total angular momentum 
L of magnetized fluid are conserved, and the angular momentum flux per unit area carried by the 
winds towards the black hole is given by —r]LB p , which is a consequence of exerting the magnetic 
torque on the disk plasma. (Now L is positive, while r] is negative for ingoing winds.) Then, the 
characteristic timescale t\ to extract the plasma angular momentum in the disk would be 

h ~ hp d r 2 d n K /(- V LB p ) ~ (r d /h)/p e n K , (47) 

where p d is the mass density in the disk with the thickness h, and (3 e is the ratio of the external 
magnetic stress BtB p at the disk surface to the local pressure in the disk. One may also consider 
the internal torque transporting angular momentum radially outward through the disk, for which 
the accretion timescale £2 is estimated to be 



*2 ~ {rd/hf/^K , 



(48) 
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and the ratio $ of the internal stress to the local pressure is usually supposed to be less than unity 
(e.g., (3i ~ 0.1). In our modeling the external torque responsible for carrying angular momentum 
into the winds should be dominant if compared with the internal torque, and we require the 
condition f3 e 3> Pih/rd, which becomes valid for h <C r^, unless the turbulent internal disk field is 
supposed to be much stronger than any ordered vertical field at the disk surface. Hence, for (3 e « 
the duration of the closed loop structure could be quite longer than the dynamical timescale fi^ 1 
particularly in a thin disk. 

In significant long term accretion of plasma beyond the timescale t\ one must consider the 
change of the closed field-line topology, unless dynamo action in the inner disk region can efficiently 
work to recover the closed loops. The process of disruption of the closed loops would be quasi- 
stationary, and our models could give a qualitative picture of the structure change as follows: 
Initially strong disk currents can sustain nearly circular loops of field lines shown in Figure 1 (i.e., 
Bq S> \B u \). As the disk currents become weaker (i.e., as the ratio Bq/\B u \ decreases) by virtue 
of the radial accretion, the loops begin to collapse towards the black hole (especially in vertical 
direction), as was shown in Figure 3. Finally, only the loops disconnected from the horizon remains 
near the inner edge (see Fig. 4 for Bq < \B U \). The loop field completely annihilates as a set of 
O-points, if the radial accretion from the inner edge continues still more (Punsly & Coroniti 1990b). 
During this process, a part of electromagnetic energy stored in loop field lines could be outward 
transferred as a Poynting flux to infinity. We expect that the magnetic field (41) to become a 
convenient initial state for analytical and numerical MHD calculations to study the problem of 
non-steady energy release in plasma accretion toward the black hole (see, e.g., Koide et al. (2000) 
for numerical investigations, in which the initial magnetic field is chosen to be the Wald solution 
in Kerr background). 

In the slowly rotating magnetospheres the closed field line plays the role of a path for the 
energy and angular momentum transfer from the disk to the black hole. [ Li (2000c) claimed that 
the transfer from the black hole to the disk occurs only for a/M > 0.36.] However, the spin-up of 
the black hole should stop when the structure shown in Figure 4 has been established. Then, the 
energy extraction from the rotating black hole becomes possible as a mechanism of energy release. 
The angular velocity Q,p of the open field lines connecting the black hole with the remote load 
may satisfy the condition £lp ~ Qh/2, which allows the twisted field with the negative toroidal 
component Bt (i-e., rj < and L < 0) to bring negative angular momentum and energy to the 
black hole. The efficiency of this Blandford-Znajek mechanism is still controversial, because the 
validity of the condition Qp « f2#/2 is based on the assumption that a causal contact between 
the black hole and the remote plasma holds. In plasma accretion onto the black hole the poloidal 
inward velocity exceeds the fast magnetosonic speed near the horizon, and no information could 
propagate outward across some critical surface (Punsly & Coroniti 1990a). In fact, the black hole- 
driven winds satisfying <C Qh have been found to exist (Punsly 1998). For such winds the main 
mechanism of energy extraction from the black hole could be the screw instability, which makes 
the twisted magnetic tubes unwind and causes a sudden transport of the magnetic free-energy into 
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the plasma's kinetic energy. For the cylindrical magnetic field shown in Figure 4 we can use the 
instability condition known as the Kruskal-Shafranov criterion (Kadomtsev 1966) of the form 

\B T \/2n5 > \B P \/S , (49) 

where 5 and S are the width and the length of the cylindrical magnetic tube. This has been proved 
by Gruzinov (1999) in the force- free magnetosphere, but it should be noted that the analysis has 
been limited to the case such that flp = and S » M (i.e., the approximation of no rotation 
and no gravity). Then, a naive application (Li 2000b) of the Kruskal-Shafranov criterion to the 
magnetosphere where Bp is roughly estimated to be r sin 6£lpB p seems to be very dubious. Here, the 
toroidal field is estimated by equation (46), which remains valid even for £lp = 0, and considering 
the region well away from the event horizon (i.e., a ~ 1), we rewrite the instability condition into 
the form 

S > \B P \5 2 /(2\ V L\) , (50) 

of which the right-hand side is constant if the magnetic flux |-B p |<5 2 in the cylinder is conserved. 
This means that if a larger angular momentum flux \r]LB p \5 2 is carried away from the black hole 
by magnetically dominated winds in the cylinder with the width S ~ M, the twisted field lines 
become screw-unstable in a shorter scale S. After the disruption of the screw-unstable structure of 
magnetic field, the rotating black hole would begin to twist again the magnetic field lines, and the 
energy release from the black hole could quasi-periodically occur. 

As was first proposed by Gruzinov (1999), the closed field lines connecting the black hole with 
the disk may be also screw-unstable. However, for the field lines closed in the disk the toroidal 
component will be mainly generated by disk rotation with the angular velocity Qk ~ Qf- Then, 
we must discuss the instability condition by including the effect of Sip, which is a problem studied 
in future works (Matsuoka, Tomimatsu, & Takahashi 2001). 

Finally it should be emphasized that the closed field-line topology presented in this paper is 
a very transient structure in rapidly rotating plasma-filled magnetospheres. The angular velocity 
ftp — &k of closed loops connecting the black hole with the disk will be smaller than the hole's 
angular velocity f2#, and the energy and angular momentum transfer from the black hole to the 
disk could spin up the disk plasma (Li 2000c). Then, as a consequence of this energy extraction 
from the black hole, outward flows of plasma from the disk would be generated, and owing to 
the frozen-in condition the field lines threading the black hole are eventually disconnected from 
the disk. The closed field lines anchored in plasma orbiting in the ergosphere would also open 
up to infinity by the strong dynamo action, as was conjectured to occur in field-aligned pulsars 
(Punsly & Coroniti 1990b). Further, the timescale ti for radial accretion due to internal torque 
may become shorter than the timescales to open up the loops closed in the disk. Then, the closed 
loops accrete toward the black hole and annihilate owing to the infinite redshift effect (Punsly & 
Coroniti 1990b). These are astrophysically interesting processes of rapid energy release, and we 
could observe transient flaring states of black hole magnetospheres, if the repeated formation of 
closed magnetic loops due to strong disk currents in the inner region occurs. 
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Fig. 1. — Poloidal magnetic field lines for the axisymmetric model of s = 2 and wq = 10. The value 
of ^ /AttM^Bq is denoted on each field line, and the structure is assumed to have the equatorial 
symmetry. Further, the central black hole with the radius r = 2M and the surrounding disk with 
the inner edge at r = 11M are displayed. We can see the magnetic connection of the black hole 
with a very narrow region in the disk. 
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Fig. 2. — Same as Fig. 1 for large scale of r, to show the asymptotic shape of field lines. 
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Fig. 3. — Poloidal magnetic field lines for the model involving a weak uniform field, i.e., 
B u = — 0.9-£>o- The magnetic field generated by the disk current is same as Fig. 1. The mag- 
netic connection between the black hole and the disk still remains, and the collimated field lines 
appear in large r. 
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Fig. 4. — Same as Fig. 3, except for the existence of a strong uniform field, i.e., B u = — l.li?o- 
Only the open field lines can thread the black hole without any connection with the disk, though 
the loop field lines still exist near the inner edge. 



